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Abstract 

The  mathematical  theory  of  propagation  for  etratified  aedia 
above  a  spherical  earth  with  an  arbitrary  scalar  radial  rariation  of  the 
index  of  refraction  is  presented  so  sg  to  be  applicable  to  both  tropo- 
Bpheric  and  ionospheric  problems.   The  method  of  traiisforni£.tion  of  con- 
tours Tised  by  Watson,  Rydbeck  and  others  for  special  variations  l6  applied 
to  the  general  case.   Sone  loathematical  questions  concerning  the  continuo\is 
spectrum,  which  are  usrially  neglected  are  examined  in  detail.  Thn  methods 
of  Furry,  Peckeris  and  others  for  the  approximate  theory  of  tropospherlc 
propagation,  which  assumes  a  flat  earth  and  a  modified  index  of  refraction, 
is  shown  to  follow  from  the  general  theory  by  ma>lng  suitable  ass^'iraptions . 
The  practical  problem  Is  reduced  to  that  of  finding  the  el^enralups  of  an 
ordinary  differential  eqiaation.   The  Vf.K.3.  method  is  applied  to  obtain 
explicit  results  for  a  ductless  atmosphere.   Tor  tlie  troposphere  a  new 
method  of  obtaining  en  approximation  in  the  case  of  a  d^jict  is  obtained 
which  appears  to  be  sulta'''le  for  quAntitative  work. 
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1,   latrodnetion. 

The  protlam  of  wtplalnln^  the  propagation  of  radio  uares   arotind  the  sur- 
face of  the  earth  has  long  been  of  great  interest  to  both  physiciste  and  nathamat- 
Icians  becsoise  of  Its  practical  importance  and  because  of  its  considerable  mathe- 
matical difficulties. 

There  are  two  fundamental  approaches  to  the  theory  of  such  propagation. 
The  first  HTrproach,  which  is  that  of  Watson  \ji}  ,   Van  der  Pol  and  Bremmer  ^2],  and  Ryd- 
beok[3j,  and  which  has  been  applied  mostly  to  lonosr^heric  problems,  assumes  the 
earth  to  be  a  perfect  sphere  with  finite  or  infinite  conductivity,  and  expresses 
the  field  due  to  a  radiating  dlpole  as  an  infinite  series  inrolvlng  spharlcpl 
Bessel  functions.   Since  the  series  converges  much  too  slowly  for  practical  purposes 
it  is  transformed  into  a  contour  Integral  which  can  be  evaluated  by  the  method  of 
residues.   To  find  the  residues  of  this  integral,  highly  complicated  asymptotic  ez- 
paaslons  for  the  spherical  Bessel  functions  must  be  used.   These  asymptotic  exoan- 
sions,  which  are  derived  by  the  method  of  stationary  phase.  Involve  the  Hankel  func- 
tions of  order  I/3  or  the  Airy  integral,  both  of  which  occur  frequently  in  diffrac- 
tion theory. 

The  second  approach  to  the  problem  is  that  used  by  Prj-ce^Uj  and  Booker  [^5]  in 
Xngland,  and  by  Turry  [6J  and  Pekeris  ["yj  In  the  United  States  while  studying  the  effect 
of  the  troposphere  on  radio  wave  propagation.   In  this  approach,  the  earth  is  assumed 
to  be  flat  but,  following  a  sugi^estion  of  Schelling  euid  Burrows  ^s],  the  index  of  re- 
fraction is  modified  so  a«  to  take  into  account  the  curvature  of  the  earth.   The  field 
due  to  a  radiating  dipole  is  now  expressed  as  an  infinite  integrad  involving  ordinary 
Bessel  functions.   This  integral,  also,  can  be  evaluated  by  the  method  of  residues 
but  it  was  realized  that  the  residues  are  really  the  successive  terns  in  an  eigen- 
function  expansion  for  the  field  of  the  dlpole.   The  propagation  problem  is  therefore 
reduced  to  the  study  of  an  eigenvalue  problem  for  an  ordinary  differential  equation. 
In  ceise  the  atmosphere  above  the  earth  is  assumed  to  have  a  constant  refractive  index, 
the  eigenfunction  expansion  is  in  terms  of  Hankel  functions  of  order  I/3,  exactly  the 
same  functloxis  which  aiipear  in  the  first  approach. 

In  this  paper  we  start  with  Maxwell's  eouationf«  and  then  present  a  complete 
exposition  of  the  theory  of  propag'  tion  in  an  inhooogeneous  atmosphere  where  the 
dielectric  constant  is  stratified  in  a  radial  direction.   The  theory  can  be  applied 
to  both  ionospheric  and  tropoipheric  problems.   The  methods  we  have  used  and  the 
results  w*  have  obtained  are  eiailar  to  those  contained  in  the  papers  by  Watson  [ij 
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Yan  d»r  Pol  and  Brenasr  [2] ,   and  In  the  books   tiy  BreinB«r(li]  and  ^dbeck   f"}].      Th« 
novelties   of  our  treataent  and  resxilta   are   the   follovin^: 

1)  The  dielectric  constant   is   considerud  to  he  a  single  function  of  posi- 
tion instead  of  different   fonctions   iA  different    regions   of   spnce,      A«   a  result,    the 
fonaulas   for  the  Herti  potential   and  the  electro-magnetic  field  hsve  a  simpler  appear- 
ance and  are  not   as   conplicated  and  cluasy  as   in  the  usual   treatment.     All   the  neces- 
sary  tranefonaations   and  approximations  can  he  carried  out   on  these   formulas   and  he- 
caU3e   of  their  siopler  appearance  it   is  much  easier  to  acquire  an  insight   into  their 
meaning.      Of  course,  when  quantitative   results   are  wanted,    it   is  necessary  to  con- 
sider the  dielectric   constant   as   composed  of  different   functions  in  different   regions, 
solve   the  expropriate   differential  equation  in  each  region  and  then  match   the   solu- 
tions   across    the  "boundary  to   ohtaln  the  complete   solution. 

2)  The   expression  of  the  Hertz  potential  as  a  contour  integral  and  the   trans- 
fonnation  of  the  contour   Integral   into  a  siun  of   residaes   eire  carefully   discussed.     There 
are  two  questions   in  this  procedure  which  are  usually  ignored.      One   is   the  question  of 
whether  the  contour  integral  over  an  Infinite  semi-circle  may   he  neglected.      Wat8on[lj 
has   shown  that   in  the   case  of  a  coixatant  atmosphere  these  integrals  vanish.      However, 
his  proof  depends  upon  the  use   of  asymptotic   fonnileis  for  the  Bessel   functions.      We 
show  in  the  Appenxlix  that,   even  for  a  completely  arhitrary  atmosphere,    these   Integrals 
vanish.      The  second  Question  is   this:      In  general,    the  contour  integral   can  be   ex- 
preeaed  as  a  sum  of  residues  plus  an  integral   around  the  imagixxary  axis;    when  can 

this   integral  be  neglected?     Vatson  [,lj  shows   th'^t,    in  the  case   of  the   constant   atmos- 
phere,   this  integral   is  negligible   compared  to   the  residue  series.     Later  writers  Qn 
this   topic  either  drop  the  integral  without  coament^or   change   the  boundary  condition 
at  the  center   of  the    earth[i^,or  express   the  Hertz  potential   in  terms   of  incident  and 
reflected  waves [2^,  which  result   only   in  a  residue   series  and  not  in  an  integral.      This 
latter  procedure  can  be  applied  to   the  general   case  but   it  is   only  in  the  case   of  the 
constant  atmosphere  that   the  integrail   over   the  imaginary  axle   does   not  appear.      In 
the  appendix  we  consider   an  arbitrary   atmosphere   and   show,    first,    th»t    if   the  earth   Is 
assumed  to   be  a  perfect   conductor,    the   integral   over   the  imaginary   axis  vanishes 
identically;    and,    second,    that   if   the   earth   is   not  aesximed  to   be   a  perfect   conductor 
but   if  its    dielectric   constant   is  assumed  to  be   large   compared  to   thst   of  the   atmosphere, 
which  is  usually   the  case   in  practice,    then   the   integral   over  the   imaginary  axis    is    small 
compared  to    the   first   term  of  the   residue   series. 
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3)     The  general   nature  of  our  study   e7nphr>slie8    the   fact   tkst    the  propaga- 
tion problan  reduces   to  the  eolutlon  of  the  eigenvalue  problem  for  an  ordinary   differ- 
ential  equation.      This   fact,   while   It   Is  clearly  recognized  In   the   flat   e^^rth  theory 
and  while  It   Is   Implicitly  contained  In  most  prerlous  work,   has  heretofore  not   been 
explicitly   stated  for   a  epherlcal   earth.      Once   this    fpct   Is   known,    It   Is   clear   that 
the  ^f.I.B.  method  should  be  applied.      By  means    of   this  method,    the    "phase    Integral" 
method   of  Bcker8ley[lOj  can  be  xinderstood  and  made  more   exact, 

k)     It    Is   shown  that   the   flat   earth  theory  follows   from  the   eixact   spherical 
earth  theory  by  making  suitable  approximations   In  the  ordln^iry  differential  equation 
mentioned  In  3),      ?rom  the  nature   of  the  approxlmatlona.    It   le   concluded  that    the 
eigenvalues  will   be  closely  approximated  but   that   the   express lonB   for  the   field   on 
the   flat   earth  theory   will   be   inaccurate  at  large   distances   above  the   earth.      These 
conclusions  agree  with  those  of  Pekerisj^Tl  in  his  analysis   of   the  accuracy  of  the 
flat   earth  approximation, 

5)     In  the  case   of  tropospherlc  propagation  the  W,K,B,   method  gives   explicit 
res^olts,      Ilrst,    if   the  atmosphere   is   inhomogeneous  but   does   not  h«ve  a  duct,    the 
field  can  be   obtained  by  a  chajage   of  scale  from  that  produced  by  a  constant   atmos- 
phere.     This   resijilt   generalizes   a  similar  resttlt  obtained  by  Bckersleyfio)  and 
Br«mner^lJ  for  a  special   fonn  of  rariatlon  of   the  dielectric   constant.      Second, 
if  the  atmosphere  does  have  a  duct,    the  differential  equation  that  must   be   considered 
will  have  two   transition  points   and  its   solution  will  depend  on  confluent  hypergeo- 
metric   functions.      We   find  that   the   differential   eouatioa  can  be  approximated  by 
one  which  has   a  double  transition  point  and  that   then  the  solution  can  be   expressed 
in  terms   of  Hankel   functions  of   order  one-fourth.      This   approach  seems   to   be   com- 
pletely new  and  suggests  a  simple,    direct   quantitative  method  for  approach  t-    the 
problem  of  propagation  in  a  duet. 

The  problem  of  propagation  in  an  Inhomogeneous   atmosphere  is   formulated 
mathematically   In  section  2  and  It  is   shown  that    the  problem  requires   the  solution 
of  a  scalar  partial   differential   equation.      In  Section  3  this    equation  is    solved 
by  the   classical  method  of  separation  of  variables   and   expansion  in  terms   of  Legendre 
polynomials.      In  Section  U  this   solution  Is   represented  as   a  contour  integral   and 
evaluated  by  the  calculus   of  residues.     These  residues   are  interpreted  in  Section  3 
as  an  expansion  of  the  solution  in  texms   of  radial   elgenfunctiens.      In  Section  6 
we  use   the  fact   that    the  wavelength  of  the  radiation  is  very   small   compared  to   the 
radius   of  the   earth  to  obtain  slimier  boundary  conditions   at   the  earth's   surface 
and  a  simpler  expression  for  the  solution.     Ve  return  in  Sectioa  7  to  the  discussion 
of  the  radial  elgenfunctlosand  we   show  how,   by  the  use  of  Langar*s   extension  of   the 


V.K.B.  method,  useful  approximations  to  the  eigenfunctlona  can  be  found, 

Th«  case  of  a  uniform  atmosphere  is  discusBed  "biy  this  method  in  Section 
8  and  the  well  knovn  resiilts  are  obtained.   Section  9  shows  that  the  case  of  the  non- 
uniform atmosphere  without  ducts  can  be  treated  almost  as  easily  as  the  case  of  a 
uniform  atmosphere.   In  Section  10  the  'flat  earth"  theory  and  the  modified  index 
of  refraction  are  shown  to  be  reasonable  approximations  to  the  exact  theory.   Sec- 
tion 11  discusses  the  case  of  an  atmosphere  with  ducts;  Section  12  sumnariies  the 
work  for  the  electric  dipole.   The  detailed  mathematical  questions  concerning  the 
countour  integral  representation  are  discusead  in  the  appendices. 
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2,      ?or7iailatioxi  of   the  Problan. 

JLgffuae   the  earth   Is   a  sphere,    radlue   a,   with  complex  dielectric   constant  € 
Introduce   spherical   polar  coordinates,    r,    e,    (J,    where  r   is   the   distance   fron    the   center 
of   the   earth,    6  is    the   co-latitude   ajid  jj  the   longitude.      We  assume   that   tho   ea'-th   ia 
imuersed  In  a  ■edium  whose    iielectric  constant,   e    ,    is    a  f'inctlon  of  r  il-ne.      The 
ornp.V-st  V- .:      r    ■■    •»    -!»»  ♦si*    ■>?>    f    'A\r.^    t    -•    f    »»    •«    rj-'jdi^^^d  'H.  ;nl-     i\\,    It   w'^ich    it 


radiating      t  a   constant   angular    frequency  •-*>  ana    Is    Irrated  at    a  dlBtunce,    1-,    from 
the  center  of   the    earth.      By   a  suitable  rotation  of  the   coordinate   syetem    .e   can 
always   locate    the   dipole   on  the  polar  eucis   as    indicated   in  the   abore  diagram. 

The  electromagnetic   field  produced  by  the  dipole  should  sp.tisfy  Maxwell's 
equations,   both  Inside   and  outside    the  earth,    and  sho'ild  possess    at    the  point    D  a 
singularity  corresponding  to   that  of  a  dipole  while   on  the   earth's   surface   the   tan- 
gential coB^onents   of   the  electric   field  S  and  the  ma^etic   field  H  must  be   con- 
tinoous. 

Let  £e  represent  the  rector   electric   field  intensity   and  He  the 

Tector  magnetic   field  Intensity  produced  by  the  dipole,    then  Maxwell's   equations 
in  M.K.S.   unite   become   the  following: 

y  X.  S  =  iu;/c  H  ,  ^  -H  "  0 

Here,   6   =6(r),    is   a   scalar   function  which  represents    the  rarlation  of  the   dielectric 
constant  of   the  medium.      Jot  r  <a,  6  (r)   is  assumed  to   be  a  constant,    ^.,  representing 
the  coi^lax  dielectric  constant  of  the  earth. 
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Vfe   shall    derive   tfte  well   knova  fact^gjthat    In  spherical    coordinates   the 
electromagnetic   field  given  ty    (1 )   can  be  obtained  from  a  Hertz  rector  wi.lch  has 
only  Its  radial  component  non-zero,    even  when  the   dielectric   constant   is   p.  function 
of  r.      Since  f  S  is   a  vector  whose  divergence  is    zero,    we  may  write 

(2)  fei  =  I  u^/tVx£   C      . 

Substituting  this   expression  for  &  in  the  curl  H   eqxiatlon,   we  get  i 

Vx  (H  -iuy«.  ^  c ;  «  0  , 

eo  th^t  the  expression  In  pprenthesis  must  be  the  gradient  of  some  scalar,  which  ve  call^. 
We  hnv.2  then. 

(3)  E  »  cu^^  C  ♦  V^  .  j 

Using  (2)  and  (3)  In  the  remaining  equations  of  (l),  we  obtain  the  following 
equations  for  C_  and  r  : 

i 
(5)  ^.Ca;^^6  c)  +  V^V>  -  0   .      .         • 

Hote   that   if  we   take   the   divergence   of  equation   (U),   we   obtain  equation   (5)   so  thst 

it   la    sufficient  to   satisfy  equation   (U), 

Let  us  express    these   eqxiatlons    in  spherical  coordinates.      Since  fr  is   a  function  j 

of  r   only   and   since   the  dlpole   Is   on  the  polar   axis,    the   electroma^nrtt ic   field  and  j 

therefore  also   the  Hertz  vector  C_  will   not    depend  on  the  longitude  0  but   only  on  I 

the  coordinates  r,    0,      We  assume  that    the  vector  £  ha.a  a  radial  component   of   the 

form  rU(r,9)  where  UCr,0)   is   some  scalar  function,   and  that   its   other  components  are  , 

*  ,  i 

zero,      Now,    the  O-component   of   KU)  iaplies   that  I 

I 


■o  that  we  may  take 
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and  then  the  r-cooponeat    of  (U)    rednoes   exactly   to   the   following: 

(6)  ZiU  *  k^  U  =•  0 

where 

and  /^  Is    the  Laplaclan  in  spherical    coordinates. 

"The   field  coioponents   can  be  obtained  from   equations   (?)   and  (3).      We   find  tJiat 

c)  r 
(8) 

The  equatlotts  In  CS)  represent  the  field  of  a  magnetic  dipol*.   There  are  similar 
expressions  for  the  field  due  to  an  electric  dipole.   The  modifications  of  the 
thiory  for  this  case  are  detailed  in  Section  12. 

Th«  problem  of  the  propagation  of  radio  waves  around  the  earth  now  has  been 
reduced  to  that  of  solving  the  scalar  wave  equation  (6).   We  know  that 

2,  ,  2 

k  vr;  =  constant  =  uJ^,>u.    ,r<a 

k^Cr)  -  UJ^6(r)yU.  .         ,  r  >  a 

and  that,  at  r  =  a,  U  and  -^s—  must  be  continuous  because  the  tangential  components 
of  the  field  given  by  (8)  must  be  continuous.   Wo  know  also  that  U  nust  be  regular 
everywhere  except  at  the  dipole,  that  is,  at  r  =  b,  0  =  0,  Tinally,  we  know  that 
the  dipole  is  sending  out  radiation  so  that  all  the  energy  is  going  to  infinity. 
This  Is  expressed  mathematically  by  the  Sommerfeld  radiation  condition.   Before 
writing  it  down,  we  remark  that  since  the  mediujn  must  eventually  became  free  space, 
the  value  of  6  will  approach  that  of  free  space  which  wo  designate  by  6  .   The 
Sommerfeld  radiation  condition  may  now  be  expressed  as  follows: 

(9) 


r->oo 


where 

(10)  k^  -u;^«s  /c 

o  0 

ik  r  / 
This   condition  means    that,   at   infinity,  U  behaves   like  •     ^  /T  »o  that  we   are 

dealing  with  an  outgoing  wacve. 
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3.      The  Series   Solution  for  UCr.e). 

The   f^ct   that  U(r,9)  possesses  a  dlpole   8in^U.<irity  at   r  »  b,   0  »  0  while 
otherwise   It    satisfies   the  wave   equation  (6)   can  he   expressed  rery  conclRely  "by 
using  the  Dir<^c  ^-function.      We  writ* 

(11).  Ziu*k2(r)U--   _/^r  -  b)  d^Sl 

2n  r"  6ln  © 

2 

where  the  denominator  2n  r  sin  9  has  been  Introduced  to  normHllze  the  solution. 

It  is  essentially  the  Jacobian  of  the  transformation  from  rectangular  to  spherical 
coordinates. 

Since  U(r,9)  ia  a  regular  function  of  (r,0)  except  at  r  =  b,  9  =  0,  we  maj" 
assume 

00 

(12) 


L^         2       n     n 


n=o 


Froffl   {\?.)  by  usin^   the  well  known  orthogonality  properties   of  the  Lefjandre  polynomials 
we   find  that 

f" 

(13)  u  (r)  =    I         U(r,e)  P„(cos  6)   sin  9  d  6 

How,  'lultiply  (11)  by  P  (cos  O)  sin  9  d  9  and  integrate  between  0  and  "n. 
The  equation  transforms  into  the  following  ordinary  differential  equation: 

1  '^^^^n^  ,  [P,s       n(n  ^  l)  1  ..  ^      /(r  -  b) 


(lU) 


or 


[1,2,  N   n(n  ■*•  1)  1 
ic  (r) T— J^ 


^         aJ-  I  2    I  n        ^  «  2 

dr        U         p    J  2  TT  r 


dr^. 


[^^(^^-  ^^]  <«.'■-  -fe^    . 


Let  us  iBTeetlgate  the  meaning  of  (lU),   If  we  integrate  it  first  fr 


oza 


r  =  r     to  r  =  b  -  h   and  than  frco  r"r^tor"b*h   ,ha  small  quantity,   we 

0          dCru^;  °  1 

find  that  — -^ has  a  Jump  of  magnitude  -     rrr*    *•  *"  g^es   from  b  -h    to  b  +  h    , 
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whlle   m     remaine   continuous.      Hote   that  the  point   r  =   a     oemses  no   trouble  'because, 

thera,    Just   as   for   every  value   of  r^lj,    the   solution  ru     and  Its   derivative  must 

n 

remain  continuous.      (These  conditions   nre,    of  cour-^e,   the   sane  as   those  which  have 

alT-eady  been  imposed  on  U(r,9).  )     The  other   conditlone   on  U,    that   It  be   rsgialar   at 

r  «  0  and  that   it   satisfy   the  Sonanerfeld  radiation  condition  (9),   imply  that  u     Is 

11 
regular  at  r  »  0  and  that 


/ 1  c-  \  i   a^  ru 


d(ru   ) 


r-^00   « 

In  order  to  scire  (lU)  we  mast  know  the  solutions  of  the  homogeneous  equation: 

dr       ^         r 

We  shall   study  this   equation,   not   only  for  integral  values   of  n  which  Rre   the  only 

ones  needed  to   solve    (14),  hut   also   for  conplex  values   of  n  which  \'ill  toe   needed 

2  2 

Inter.      In  case  k  (r)   is   a  constant,   k     say,    the   solutions   of   (16)   are  arbitrary 

linear  combinations    of  r  '      J  ^,  /^(k  r)  and  r  '      H\.,  /^  (k  r).        The  function 
,,  nTl/2     0  a*-l/2       o 

r"^'      J  ^i/p^^Q^^   1-s   regular   at   the  origin  If  real  part  of  n  >  -  1  while 

1/?      (1 )  ? 

the  f\xnctioa       r  '      H  ^i/J'^^^   satisfies   the   radiation  condition.      If  k     Is  not  a 

constant,   the   solutions   of  (I6)  will   behatve  appro xijnately  like   the  Bessel    functions 

2 
but   their  exact  behavior  will   depend  on  the  way  k     varies  with  r. 

We   shall    Introduce  three  particular  solutions   of  (I6).      I'et   rr   (r)  be  that 

(1)  ° 

solution  vrtilch  Is   regular  at  r  =  0,   while  rw  (r)   is    that    solution  which  fox  large 

r    is  aeymptotlcally  e^'   ,   and  rv^     (r)  is  that  solution  which  for  large  r  Is 

ftsyaptotlcally  e"^  .     Bote  that  the  fanction  rw       (r)   satisfies   the  radiation  coa- 

n 

ditloa  (15).   Tor  slBpliclty  of  writing,  we  shall  write  It  as  rw  (r)  when  there  !• 

n 

no  possibility  of  eonfusloa. 

Since  the  differential  equation  (I6)  is  unchnnged  when  n  is  replaced  by 
-n  -  1,  we  see  that 


>rlgln  will  in^jly  that  ▼  (r)  behaves  like  r'*(n  >  0)  at  the  origin  and  this  con- 


This  relation  is  not  true  for  ▼  (r)  because  the  condition  of  regularity  nt  the 

,TJ    oennves  i.'Licm    r  ' 

ditlon  Is  changed  when  n  Is   replaced  by  -  n  -   1.      It   sliould  also  be  noticed  that 

2  2 

for  k     =  k     ,   a  constant,  while 
0 
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on  th»   other  hand,    because  of  the  behprlor   at   Infinity,    we  hare 


and 


n 


.^.-WU    .-i(^l)n/a(,^,,i/^  ,(.)^^  ,,^^j      _ 


We   shall   now  prove   the  well-knovm  fact   that   the  Uronskian  of  any  two   solutions 
of  (lb)   is   a  constant   independent   of  r.     Consider,    for  example,    the  two  oolutlona 


rr  (r)  and  rw  (r),   We  hare 
n         n 


^.  .   ai^] 


^(-„'*h'-    ^^[("^J'O 


;rw.)  .  0 
n 


^  (,.„).  [.^.  sis-)], 

dr         ^  T        ■* 

Hultiply  the  first  equation  by  i^.,.  the  second  by  rv  and  subtract.   We  get 

« 

(rw  )  -^  (rv^)  -  (rr^)  1^  (rw  )  =  0 
°  dr^    ""      ^  dr^    "^ 

Integrate  this  equation  from  r  =  r,  to  r  "  r_,  emd  we  hare 

SO  that  the  Wronskian 

(17)  W^  =  rw  (r)  1-  (rr^)  -  rr  (r)  A.  (rw  ) 

is  independent  of  the  value  of  r  but  depends  only  on  the  value  of  n. 

Now  that  we  have  soae  knowledge  of  the  solutions  of  (16),  we  can  prove  that 
the  aolution  of  {Ik)   is  given  by  the  following  formlaa: 


r-r^ 
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a 


(IS) 


rv^(bV(r) 
n 


r  >  b 


First,  it  is  clear  that  both  expreseioa*  in  (18)  satisfy  (lU)  for  r^b.   Second,  at 

r  =  b,  tha  function  ru  is  continuous  and  its  ieriratiTe  has  a  jump  of  magnltud* 

-l/2nb.   Third,  from  the  definition  of  the  function*  t  (r)  and  w  (r).  the  function 
'  n        n 

u  is  regular  at  r  =  0  and  satisfies  the  radiation  condition  (15),   This  shows  that 
n 

ru  (r)  really  is  a  solution  of  (14)  satisfying  all  the  other  conditions. 

Finally,  we  can  obtain  a  representation  for  U(r,e)  fro«  equation  (12).   We 

have 

00   T  (r)w  (b)   -  ., 
U(r,e)  -  ^   y-      '^y  °     ^  P^  (cos  9)   .  r  <  b 


(19) 


1_  oo_  V^)w^(r) 
3"  Z_-     W, 


~ip^  (cos  e)    ,  r  >  b 
2    n 


This  is  the  solution  of  the  problem  in  its  most  (general  form  for  an  arbitrary  variation 
of  the  dielectric  constant.   Special  cases  of  this  are  well  known  for  the  cases  in  which 
the  ordinary  differential  equation  (lb)  can  be  solved  in  terms  of  known  functions;  for 
exairple,  BydbeckjjJ  has  treated  the  pro'olea  where  ^  (r)  is  atsuaed  to  VHiy  parabolically. 
It  is  to  "be  noted  that  thia  solution  (I9)  can  be  used  for  conputation  only  if  ka  is 
small  coapared  to  unit/,   In  the  following  sections  this  solution  will  be  tranaforned 
80  that  it  can  be  readily  computed  for  large  values  of  ka. 

It  is  iaterestizig  to  consider  some  special  cases  of  (I9).   Suppose,  first,  that 

2    2 

the  mddiuB  is  uniform  and  that  there  is  no  earth  so  that  k  «  k  .  a  constant.  t-i«n 

oL  • 

using  the  standard  notation  for  the  spherical  Bessel  ftinction  ri3J 


-IP- 


^n^^^  "  V^'^ 


(r)=l'^*^h^lN).r) 


»n 


f onula  C19J  bacoBft* 

Ik  H  ^  ^ 

0 

(20; 

'^     ^    YL     ^"^"'^   ^^^V^  ^1'^    ^^'^  ^n^'°^  e)    .      r  >  b  . 

0 

the  Classical  formula  for  the  dlpole  [ l3j  t  where  R  is  the  distance  fraa  the  point  with 

coordinates  (r.O)  to  the  source  point  ('b,0).   Suppose  now  that  the  earth  is  present 

2    2         2    2 
so  that  for  r  <  a,  k  =  k.,  r  >  a,  k  «=  k  ;  then 

n  nno  nn  o 

w   (r)  =  i'^"^  h^^^    (k  r)  ,  r  >  a 

n  no' 

where  ^     ,    0       are   to  \)e   datemined  by  the   fact    that  u  (r)  and  t"0.  (r)   are 
n         n  n  or     n 

oontimioas  at   r  =  a.      Troin  the   definition  (I7  J  of  the  Wronsklan,    it   follows   that 

w  .  i'""^  b^v   k  fh^^^  (k  b)  j'  (k  b)  -  h^^^'  (k  b)  J    (kb)l 

n  nol.no         no  n  o*'noj 


1»  X    /k 
n'    o 
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X  is   foiind  froB  the   following  equation*: 

n.     X  n     n     0  &     n  o 

k,    j'    (k,a)   =    y„   k     j'    (k  a)  -^  r      k     h^^^'    (k  a) 
j.nl                none  non  o 


We  have 


k.    j'    (k-a)  kj^^   (k  a)  -  k^   J      (k.a)  h^^  ^   (k  a) 
y      _       1     n       1 & 0  o     n       X n o 

^"  hl^^   (k  a)   j'    (k  •)  -  h^^^'    (k  a)  J^  (k  a) 

n  0         n       0  n  0         n       o 


Bc   that 


(a)  ^   „  ,nn       2  r       '  (    ^^  ^(1)  ^    ^j  _  (      )  ^(1)  ^      1 

n  o^i.a.i&  0  onin  oj 

finally. 


00  J    {k,r)lii^\k  b)i°"'S  (co8  9) 


_      1 


(2^1)       "     "       " ^ ii .      r<a 

V 
o  n 


(2?)  ,00  Y  J   (k  r>i*  h^^\k  r)     _^,      ,.  . 

^  »     ^    y     (2un)       ^'^  °     %       -^  ^ 2_  i""'  4'^(V)  P,(ccs  e).a  <  r  < 


0 


00  jr  J    U  b^J'h^^^Ck  b) 


00  AJVKg/TJu  \KO/  .-         /-\ 


n 

0 


Sotice   that 


n'    n         0         * 

80   that  'by  using  (20),    (22)   can  he  vritt«n  as  the   gram  of  two   terms: 


-lU- 


ik  S 

(23) 


*— —  n 


The  flrct  term  Ib  due  to  the  dlpole  while  the  second  term  is  due  to  the  scattering  of  the 
dipole  field  hy  the  earth. 

In  oaee  the  eprth  Is  assumed  to  te  a  perfect  conductor  "80  that  Ic,  =  00,  then 
J  (k_a)  »  0  and  from  the  ccntinulty  condition 

^n^-^n-JnV^/^i'^^V^    • 

We  find  then  that 

ik  R   (k  h) 

t7(r.e)-  ^  y"(?n*l)  n  (k  r)h^^^k  a)V^^(k  r)  J  (k  r)l   f,  ,   °   P  (cose)  a<r<t 

n    0 

*——  *-  *  h   (k  a; 

n    0 
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U.    ContottT   Intagral  Representation  for  U(r.O) 

7ollowln^  a  method  due   to  Wateon   [l]     .   we   shall  ehow   that   the   serlee    (19) 
can  'be  ezpreeeed  aa  a  contonr   Integral.      Consider   the   eerlet 

(25)  s  -  i:    2 

o 


2n+l 


a     P      (cob   9) 
n     n 


whf>re  a     Ip  easumed   to  be  an  analytic   fxinction  of  n.      We  ehall  show   that  S   cen  be 
n 

represented   by  the  following   Integral: 


(26) 


2tt1  J- 


(t+1/2)  a^  P^    (-cos   ©) 

sin   TT    t 


dt 


where  C  =   C.    +  C-   is  a  contour   that  starts  at    oo  -   1 d     in   the   t-plane,   goes  below 
the   reel  axis   to   t   =  -   l/2  and  then  above   the  real  axis   to    oo  ■♦•   1  (T    . 


/-  plane 
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Slnce  a  and  P   (-cos  0}  are  analytic  functions  of  t,  the  only  Blngnl- 

V  t 

arltles  of  the  Integrand  are  poles  at  those  Talues  of  t  inside  C  for  which  sin  -nt  =  0, 
that  It, 

t  »  0.  1,  2.  ... 

Since  the  contour  C  Is  descri'bed  In  the  clockwise  direction,  the  integral  I  is  equal 
to  the  negative  sum  of  the  residues  at  the  poles,  so  that 

oa       (n+l/2)a  P  (-cos  0)  oo 

(27)         I  -  -  21  ^  ^„;  °  -  -  ^  Z  (°"^l/2)  a  Pn  (  =  °«  ®) 

^r~  Tr  cos  TT  n  TT  *~-         n  n 


because 


P  (-C08  «)  »  (-1)°  P   (cos  0) 
n  Q 


Tormula  (2?)  shows  that 

S  =  -TT  I 

If  this  procedure  is  applied  to  (19) t  we  find  that 


,     /•  (t+l/2)T.(^)w.(r)   P.(-co8  0) 


C 


where  C  is  the  contour  specified  abore.   On  C^ ,  the  part  of  C  below  the  real  aiia, 
replace  t  by  -t-1,  then  C   is  transforoed  Into  C.  and  we  have 


(29) 


.        .  if  (Ul/2)r,^_^(b)w_^.^(r)   P,,.,(-co»  9) 

^^'^^^'-     Hii^J^  w— slnTTt ^' 

^   /*     (t+l/2)T^(b)w^(r)    P^(-coe  ©) 
"  h^Jc^  Wl  sin  TT  t     ^^        ' 


i 
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We  Bhftll  rotate  C^  and  C-  to  that  they  go  around  the  upper  half  of  the  line, 
laa^lnfiry  pert  of  t  =  -  1/2.  To  Justify  this,  it  Is  necessary  to  show  that  the  inte- 
grals   over 


Im(t)'--L 


t '  plane 


the  infinite  guart'-r-circlee  in  the  first  and  second  quadrants  go  to  lero.   ftii«  :  • 
discussed  in  the  appendix. 

When  the  contour  is  rotated  in  the  desired  inanner,  it  may  pass  across  some 
poles  of  the  integrand.   We  find  then  th«.t  U(r,©)  as  given  by  (29)  can  be  written  aa 
the  SUE  of  residues  at  the  poles  plus  two  integrals  orer  the  upper  half  of  the  line 
lB(t)  -  -  1/2,   Put  t  -  -  1/2  -►  if  eJid  use  the  fact  [l4]  that 


(30) 


P^(-coe  e)  =  P_t-1  ^"*'**"  *^   • 


we  obtain 


(31)   n(r.Q).  r  residues.  5^5- J    —f^osh  ^r 

•  vo  * 


V-l/2^^^^r-l/2^''> 


'it  -1/2 


(b)w 


^-ir-i/2^^^''-ir  -1/2 

"-IT  -1/2 


,(r) 


dr 


From  (28)  it  is  dear  that  the  poles  of  the  integrand  will  be  at  the 
points  where 

W^  -  0   . 

Call  thiB*!  points  t^ ,  t    ...  ;  then  th<*  reeldue  at  any  pole^t  ,  eay,  is 


(32) 


Here 


-IS- 


(t  >l/2)v      (b)v      (r)P      (-C03  0) 

^  ^1  ^1  ^1 

W;      Bin  TT  t. 


(33) 


w>    -    —   w 

t  dt        t 

J  't-t 


We    thuB  have   the  final  result: 


D(r,Q)   =^ 


(t  >   l/2)v      (b)w      (r)P      (-coe   C) 

.J !i__j J 


^        Bin  n  t 


'^/: 


00  '^Pi-r-_i/2^-<=<'«  ''^ 

-COEh  TIT 


.(to)w,^     wo(r) 


^r-i/2^''^V-i/2^'^^      ^-ir.i/2^^^^-ir-i/2 


(b)w   ,^   Wo(r) 


'ir  -1/2 


'-ir-i/2 


dr  . 
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5.  The  BlfcgnraluB  P^oblea. 

The  residiie  In  (32)  caa  be  aore  ueefull/  Interpreted  In  another  way.   Since 


the  Wroneklan, 


W^   -  0   . 


the  two  functions  t   (r)  and  w   (r)  are  not  independent;  that  is,  there  exist  con- 


stants  a.  such  that 


^  ^'^  '  "^i   "t  ^'^    • 


Thl«  equation  shows   that   the  function  w      (r),   "besides   satisfying   the   radifttlon 

condition  at   infinity,   also   BStltfies    the  condition   of   re^larity  pt   r  =   0.      In 

other  words,    t     is  en  eigenvalue  and  w      (r)   is   an  elgenf ijnction   of   the  following 

^  *J 

problem: 

rind  th«  values  of  t  for  which  T(r),  the  solution  of  the  following  differential 
equation: 

(3M  -^  *  P  -  i^1  (rr)  .  0    , 

dr  r    •' 

is   regular  at  r  »  0  and  satisfies   the  radiation  condition   (I5)  at   infinity,    that 
is, 

Itm     I    dCrr)  .,         1        . 

^         —J — ^     -   Ik  rr     «  0 
r->oo'      dr  0      ' 

Pttt 

(35)  P       =/      w^     (r)^     4r 

J        ^  0         J 

1/2 

so   that  w      (r)^,  i"  a  noraaliied  eigenfunction   of    (3^);      then   the   residue  given 

tj  J 

Ijy   (32)  will  be  shown   to  be  equal   to 
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(36) 


w      Cb)  w      (r)   P      (-C08   ©) 
eln  IT  t, 


2  pj  ...  ..    .J 

This  resnlt  will  be   obtnined  from  a  consideration  of   the  following  two  eqijAtlone: 


d^Crv^) 

t(t+l)  1 

t.(t  +1). 

_J — J 

dr^ 
d^nr,    ) 

]  <-t> 


=      0 


dr 


■]   (iv^    )     =0      . 


Here  t  hee  any  valiie  and  Is   not  an  eigenTalne,      Mnltiply  the  first  equation  by 
rw        ,    the   second  by  rr     subtract  and    integrate  with  respect   to  r  from  0  to    oo. 


We  obtain 


07) 


(t-t   )(t+t  -t-l) 


/oo 
0  J 


dr  =   nr^     f?  ^^t.^   "  "'t ,     d?  ^"t^ 


00 


Notice   the   term  at   «ero  drops   out  because  by  ceflnition   ▼.    and  w.      are    re^jilar  et   r 

J 

Consider  the  following  determinant: 


rw. 


TV^ 


IV. 


rv. 


J 


fe<-t^   ^-'t^   ^-t,^ 


It    is   obriouely  zero   tince   the   elements    of   the  first   two  rows  are   identicfll.      When 
it   is  expanded  in  terms   of  the  elemente   of   the  first   rov,   \rm  hare 


rv. 


-'-t  ["t^fF^-t^-^^tlF^^tjO    "    • 
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Now,  99  T   approaches  infinity,  the  firet  bracket  approaches  zero  since  rv  and 

w 

rv   both  satisfy  the  same  concHtlon  et  infinity.   This  shews  that 


d 


't.  ["t  h  ^"t^  -  ^tfe  ^^tO 
■J 


80   thfit    (37)   can  "be  written  as  follows: 

w.    (r) 


/*  *  *J  ^t 

J  ^  ^  \  *^  ~-  Tfir  (t-t  )(t*t,^i) 


As   t  approaches   t.,   w     approaches  w.      and  v     epproeches  a.  w       while  the  right 

side   of  the   above   equation  becomes   indeterminate.      If   the   right   side   is   evaluated 
by  L'Hdpital's  rule,  we  find   that 

a,         w.    (r)     dr  = 


I  \ 


i        t/  '  2t>i 


80  that 


i  (2t_,+l)aj 


and  U8in£   this   in    (32)  we  get    (36)    . 

When   (36)    Is   substituted   in    (31)  we  get  finally  the  res-ilt   that 


r^   (b)w^   (r)       Pt   ^"*^°*  *^ 


(38)  U(r.©) 


2rt  .  sin  -n  t, 


-2&- 


where  I  io  the  following  integral: 


I  = 


_^  r  00     rP^^  ,i/2(-<^°«  ^)      [Mr -1/2^^ Kir -1/2^''^ 

WlJ^  coah  nr  L  *-lr-l/2 


^r-i/2^^Nr-i/2^^)  1 


dr 


ir  -1/2 


In  the  appendix  It  will  he  shown  that,  if  the  earth  is  assumed  to  be  ] 
a  perfect  conductor,  then  I  vaaishes  identically.  If  the  earth  is  not  essiimed  to  ■ 
he  a  perfect  conductor,  I  will  not  1)6  cero.  However,  It  will  he  shown  that  I  is  ] 
small  compared  to  the  first  term  of  the  series  in  (38)  a»  long  fs  the  dielectric  ' 
constant  of  the  earth  is  large  compared  to  the  dielectric  constant  of  the  atmos-  i 
phere . 


-?■> 


b,     i   Slcpler  yprm   for   the  Herti  Potential. 

So   far,    we  hsre  not  yet  nade  use   of  the    fact    th«tthe  dielectric   conetant   of 
th«   earth  18   a  constant,    nor   of   the    fact    that,  for   the  wavelen^the   we   axe    Interested 
In.ka  Is  verv   large,    of  the   order   of  magnitude   of  10  .      We   shall  use   these   facts 
to   Binplify   the  expansion  given  in  i'^S). 

It  will  "be  seen   in  the  later   sections   that   the   zeros     of   the  Wronskian,    t., 
will   all  be  Tery  Inxge,    of  the   order  of  mrignitude   of  lea.      We    are   therefore   Justified 
in  using  an  asymptotic  formula  for  the  Legendre   functions   that   ^pear   in  (38).      Ve 
have  {[ll*] 

\    1    /  ^ 


for  large  values   of    v  ,    so   that 

^^^  slnTTt       "  sinTT   t  '^"    (nUj^l^sin  0  /  * 


2         2 

Haxt,    since  for  r  <  a,   Ic     •=  t    "  constant,   we  know-   that  the  solutions   of   (3U) 


are 


rr  =  r 


"'  •'.n/s  <^'>     -   ■•''''  <ii/e  O^i')-    '  < 


The  function  ▼^.(r)  was  defined  to  "be  that  Bclutlon  of  (3^)  which  was  regular  at 
t 

r  =  0.   It  mast,  therefore,  "be  proportional  t  ;  the  function  ^^^-t/n^^^^  ^°^   r  <■  a. 
Since  the  solution  of  (3U)  emd  its  first  derivative  must  he  continuous  at  r  =  a, 
this  requires  that 


(la)      l£^  t    I    .  ax 

rv 


Again,  the  right  side  can  be  siapllfled  "by  using  the  asynptotlc  fonaulas. 
Put 


hi'"'  ^n/.  '^-') 


g(r)  -  — YT,- 


-2U- 
to  that 

'   ''tn/2  ^''i'^  °  • 

where  C  ie  a  niitable  constant.   Substituting  this  exprestlon  In  the  differential 

1/2 
equation  for  r  '  '^*^i/p   (^*'^t  ^a*  l-^t 

ve  get 


(US)  ^    -     «=  »  icf  .     ^^f^      -  0 


-2  ^  ^2       tUn) 


2         2 

Since  t.   "U»     f -At   v/here     f,,    the  dielectric  constant   of   the   earth,   is  very  large. 


we  nay  neglect  ^     in  (H2)  and  get   th.e  approj^iaate   formula: 

g(r)  "*  1    y  k^  -     ^ 


2 

r 

To   detemine  whether  it  is  a  plus   or  ainus   sign,   we  notice   that    (:  .    and  so,    also 
t  ,   has  a  large  positiTe  Imaginary  part   so   that,    hy  the  usual    Rsymrjtotlc   formula, 

''^' Vi/2  V^^/l   -"[V-   i<tn/2)-   l\ 

-ik^r 
^Cg  e 

v^ere  C      is  independent   of  r.     Then 


siL^lViZaV) 


-?5- 


and  80  g(r),   vhlch  is  a  more  accurate  expression  for   the  ratio   of  the  derlvatlTe 
to  the   fanction  anist   have   the  mimis    sign,    that    is. 


g(r)  -  -  1^    k^  -    - 


Tvnl 


2 

r 


Ve  shall  see  later  that  t  'Vak(a).   Because  of  this. 


g(a)  =  -  ly  k^  - 


k^ 


a  quantity  independent  of  t.   Bow  (Ul)  can  be  written  as  follows: 


^:^l..  --^r<^ 


This   result  holds   for  all  ralucs   of  t,    if  t    is   approximately  ak:(a). 

Consider   the   eigenvalues  t.    for  which  the   ftuictions  t   (r)   and  w.  (r)   coincide. 
Then  (U3)  gl^es  us   the  following  eigenvalue  equation: 


*    rwt,   (r) 


We  shall  discuss  its  solution  in  a  later  section. 

There  ie  one  nore  quantity  in  the  expansion  (38)  that  will  he  simplified. 
It  is  the  quantity  0.  or  essentially  V  ,   Since  the  Vfronskian  is  constant  for  all 
values  of  r,  we  have 

^  =  ~t  h  ^'^t^       -  "^t  h  ^"'t^ 


and 


'^(rwj 


V 


i-   (rv  )  -^  rw.i-^  (rr  ) 


t  "  "TT—  d7   ^"t^  '  "'tif  dF  ^'^t 


-  d?  ^"t^  -  "t  -fr  T7  ^^t^ 


9t 
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all  eTaluated  at  r  =  a.   Now,  8inc«  the  right  side  of  (U3)  is  Indepeadent  of  t 
haTa 


we 


at  TT^ 


=  0 


r»  a 


or 


t  at  a.r       t 


Tt    dF  ^""t^j 


When  t  ■  t  ,  an  elgenralue,  the  function  t  (r)  is  Identical  with  w.  (r)  and  then 


Q      d 


"'t.  ^Tt  ^^"'t.^  - 


■^  ~  (rr.  )  =  0 
J 


Usin^  this  result  in  the  formila  for  W   we  find  that 

V 


9'~,  ) 


-'^K  -It  h  '~t  )l     • 


If  V     (a)^0,   which  V   (^)  will  always  be  the  case   if  lt^j»^oo,    that   ie,    the   eprth  is 
not  a  perfect   conductor,    then  we  maj-'  write   the  above  fonsrala  as   follov/e 


(U5) 


,2      "3        dF^'^tj^ 


^"^t/      Tt         rw, 
J  t 


(aw,    r     M 


vhere  M  ie  defined  b/  the  last  equation. 

Combining  the  results  we  hare  obtained  in  (U5),  Ckk),    (I+2)  and  (38),  we  get, 
finally,  the  desired  expression: 


-27- 
,1/2 


['(2^1)]  '        i(Yi/2)  e     "t^^^^^t^^^^ 

[v  sin  0  J  •  ' 


ii(r.e)  =  -    y~       .. ^ 


(U6) 


„      :;:1_     ^- —  1  J  £ J  J 

br     2L-.  r     8ii  0  I  N  aw     {&)     aw     (a^ 

J      *■  -^  J  J 

Here,    bw.    ("bj/aw.    (a)   la    the  height   gain  factor   for  the   transmitter  and  rw      (r)/aw     (a) 
Is    the  height  gain  factor  for  the   recelTer. 
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7,     The  Langer  Aayptotic  ypmala*. 

The  prerloTis   sections  hare   shown  that   the  problem  of  determining  the  Hertz 
rector  TJ(r,9)   and   thus  the  electrojia^etic   field   depends  upnn  the   solution  of  the 
eigenvalue  problem  for  the  differential   equation  (3^).      In  general,    this  problem 
has   to  he  sdLred  by  some  ^proxlmation.      In  this    section  we   shall  discuss  by  a  method 
due  to  Langer  [^9]    *^®  well-known  W.E.B.   approximation  to  the   Bolutiona  of  a  second- 
order  differential   equation  and  in  later  sections  we  shall  use   these  approximations 
to  obtain  the   eigenvalues   of   (3^)*  which  we  write  as  follows: 

.2 


CU7)  ^-^    -p^'>   -     •^](rv)  =  0. 


Here,  we  have  i*it  \  for  t(t*  l)  . 


Consider  a  second  order  differential  equation  of  the  form: 
dr' 


(14a)  -^  -^    [O^  f(r)  *g(r)]s  =  0 


where  0  is  assumed  to  be  a  large  constant  and  where 

f(r)-  (r-r  )•"  f  (r)   , 


o<    1b   a  constant  and  f  (r   )/0,      Put 

.  1/2 

(U9)  r-r^-h(r)  .  >A,-<r(h(r)) 

where  h(T)  Is   so  chosen  that 

(50)  h'cr)^  f  (r)  -  1 

and  T"  "  0  for  r  =  r   .     Under  the  transformations   (U9)   equation  (U8)  "becomes 

0 

^^  dl^      I  ivC-^  *2)¥  '' 

\rtxere  4^  (T)   ie   small   for  large  values   of  i)   .     The  solution!   of   (5l)   are,   approxi- 
mately,  linear  combinations  of 
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(52) 


T1/2hU)    (jd   andrV2H^2)    (;,r) 


where 


?^  =    («  ♦  2> 


-1 


We  shall  see  later  that^for  oar  purpoeea^  oaly  the  caeeE  o<  =  1  and  (2^  =  2 
are  of  intt  rest.   The  case  <<  =  1  ccrresponde  to  the  situation  where  the  atmosphere 
does  not  have  a  duct  or  a  temperature  Inversion.   The  cas*  <^  "  2  vlll  correspond 
to  the  situation  -.'here  the  ftaiosphere  has  a  duct. 

In  case  tK  *=  1 ,  we  rewrite  equation  (U7)  as  follows: 


dr 


L^«^o  I  ()cr}\ 


(rv)  -  0 


This   equation  is  of   the  saae   form  as    (U6)  if  we  make  k  a  correspond  to  ?    and  put 


^V   /  (kr; 


.       g(r)  =  0 


To  use   the  fornulas    (5?)  we  oust  put  from  (5O) 


dr         h^Tfy        L^^^<J 


or 


(53) 


=i;['^i'-^i 


1/2 


dr 


Here,  r  is  taken  to  he  a  zero  of  f(r)  so  that 
'   0 
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0  0 

Ii'  we   chango  also   the   dependent  variable   from  v(r>   to  <r  (T)  by   the  transformation 
indicated  In  (U9): 

-1/U 

r     I.  "1    1/9  l/^rd-^X"-  A/a"-     I 

rr  =     <r 


0 

then,    according  to   (^2),   we  hare 

-1/; 


i'-''' '  ^''"  WY -  ^]'' '  Of  >v« 


"'t 


In  particTilar,  the  solution  which  satisfies  the  Sommerfeld  radiation  condition  at 
infinity  will  "be  given  hy  the  formula: 

-1/U 

0 


It   ahould  be  noticed  that    the  'ilKB  method  and  the  Langer   formuleus  are  not 
TeJ-id  if  f(r)   in  {h8)   is  dlacontiimous.      Thia  meana   that   fomula  (5^)   is  not   valid 
for  r  <  a.     However,    we  have  seen  already   that   the  solution  of  (3^)  for  r  <  a 
hehaves   like  the  Bessel   function  ^ *^i f r^^^^ ^ »     Using  the  approximate  formula  (5'+) 
and  neglecting   terma   of  order  (ka)~^,   we  obtain  a   simpler  fom  for  the   eigenvalue 
equntion  (U^)«      It  becomes 

'"'        It  ■'o )  -  i.yi    ^i^r^  ^  ^,  V  M^i 

Here  T      1^   the  value   of  T  corresponding  to  r  =  a.      Hote,   that,    if  the  eeirth  is 
assumed  to  be  a  perfect  conductor  so   that  ki    =    00,    then  equation  (55)   reduces   to 
the  following  siiqsle  form: 

(56)  z[y^  (v^a^  -  ° 


-31- 


Slnce  the  rHtlo 


^/3  ''^    /    (1) 


will   appear  eo    frequently   In  our  analyslB,    we   shall  use   for  It   the   Birnuler   notation 
Z(x),      Equation  (55)   can  now  be  wrttten  as   follows: 


0-^'^ 


1/2 


Z   (k  ar  )  = 
o     a 


^  f-M 


Put 


I      Z(t)   dtl 


where  C^    and  o(  are  suitably  chosen  constants.     From  the  differential   eauation 
eatisfitid  by  y 

y"  f    -  y'  *  (1 i^)  y  »  0 

> 
we  find  that  Z  satisfies  the  Biccatl  equation: 


I    2    1  1 

2  +Z'^+  i.  Z*l„  -±- 


so  that  for  large  x 


(57) 


Z'<v-  1  -  Z^  -  -  Z 

X 


-3^ 


8,      The  Case   of  a  Uolform  Atmoaphere. 


We   shall   dlscuas   first   the    case    of  a  unifono  atmosphere    In  order   to  more    easily 
understand  the   case   of  a  non-uniform  atmosphere.      The  methods   and  res\ilt3   of   this 

section  sre  classical  juj.       Suppose   that    the   atmosohera   is  uniform  so  thnt 

2  2 

k     =  k     ■=  constant   and  also   that   the   earth  is  a  perfect  conductor.      The  solutions   of 

,(1) 


are 


Hi^^(f)»0 


y  ^  -   2.38   .^''    .  Fi  *  5.6   e-^S     f     »  g.65   e-"\      etc 


The  equation  (36)  now  hecooes  the  follovin^: 

a 
^59^  k^  /     I  1  -    -•  ^  I      dr  =  r_  -  e  "  ^n 


/    [1.  -^]'^'  dr=  E   -e^^r, 
^r/     (k^r)2^ 


where    f       Is    defined  "by   the  last   equation. 
J   n 


n 
Put 


in  (59),   Since  k  r^  =  X  '  ^  ,  it  reduces  to  this: 


o      0 

(a 


'rfhen  ve  eriraluate  the   Integral  we  get,    finally,    the  eigBinralue   eouetlon, 
(a)  (p  f  .  1)^/2  -  arctan  i^\  -  1)^^  =  ."^^    T^  X~V2     . 

In  this   equation   P       is  a  function  of    \   defined  by  (60)« 

*  1/2 

We  shall   find  that  \  '      is   of   the  order  of  magnitude  of   (k  a)   and  therefore 

the   right   side  of   (60)  will  he  rery   near  to    zero  which   Implies   that    p      will  he 

close   to   one. 

Put 


-3> 
thau  (61)  'becoaes 


J  o  a  no 


or 


(62)   •  J-.  I  (3  7-^2/3    ,-2"i/3   . 

Sach  value   of   f      determines  a  value   of    ^  ,    then  a  value  of    0    ,    and  by   (60)   a  v^lue 

Q  1/2  ^ 

of  ^>  ,      Since   t  Is   approximately    \  '    ,    v/e  get   finally  the  well-knnown  expression  fill 

for  th^   eigenvalue: 

(63)  t,-V-  I  (V>''^  (3r„)''' .™''  . 

Kote  that  (62)  shows  that  T  is  of  the  order  of  magnitude  one  and  since  (k  a)**  '  ' 

0 

is   saadl,    our  assumption  that   P       Is  close  to   one   is  .1u   tlfled.      It    is   rIsd    import- 
ant  to  notice   thax    the  integral  (bOa)  can  b«  easily  evaluated  approximately. 
Since  Q       Is  very   close   to  one,    we  ma^   proceed  as   follows; 


rea  1/2  pH         ^       1/2  ,  ,.   fPa  .  ,^ 


dp 


=     1     23/2   (  p ^-1)3/2  .     1    (2^)3/2  (j^^^j-l 

which  is  the  same  result  that  was  obtained  by  an  exact  integration. 

If  we  drop  the  assumption  that  the  earth  is  a  perfect  conductor,  then  we  must 
solve  the  eigenvalue  equ=^tion  (55^  instead  of  the  equation 


■d) 


"in'V 


0  .   , 


Using  the  same  approxlaatione  and  notations  as  before,  we  find  that 


-3V 
and  th.it 


0  o 


8o    that    (^5)  becomes    the  following: 

wnerd  »^   is    ieflned   o:'    th-   l-^st   equation.      Mathols    for  solving   this   eou^tion  are 
di->cus?ed  tiy^  5reianer   fill  , 

V.-e   shall   no-./   find  a   simple  auprozinate   formula  for   the    quantity  M  „hlch  a.-jpears 
In  (US)   and  then  w^   can  write  down  the  complete   fonnola  for  the  Eertz  potential   as 
given  -by   (Ub;.      From   ik'^),    (57)   and  (61+J   we   find  that 


D  t  rw 

tj  la 


^0  v^''/''^t^^^^'/'4^2^^'^'/3}] 


How,  \>y  (60) 

t(t*l)  =  A    =   ^V^Va"^  "  ^V^^  [l"'    J*(k^ar2/3j  -2 
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80   th'it,   approilaately 


t  =  k  all  -^  ^(k  a)-^/-^!  -^  =  k  a  -    5'(k  a) 

O     ^  0  J  0  0 


1/3 


and  thea 


.IL.^-  (kar'/3 

U3ing   tills   result   lii   the  expreealon  for  M  given  above,   we   find,    finally,    that 

(65)  M  »  k^  ^^c*'"''^^^  r^^*  ^^  ^V^^^"^  •»->j(k^a)^/3<r-^J    . 

Coni^ider  the  case   where  both   the  trar.smltter  and   receiver   are   en  the   earth's 
surface  EC   that  b  =  a  ==Lnd  r  =  a.      Substituting  (65)  and  (63)   iii  (U6)  we  get 


i 


ik  a6 

0 


Herp  T'.    is   the  root    of   (6U),     Kote   that  afi  =  d  is    the  dlstanc*  from  the  transmitter  to 
tha. receiver   ec   that 

ik  d 

u(a.e)  =  -    5^  yTSuX     f  (?^) 


wnere 

>i/3 


7^  =  (V^      ® 

J 


'   ^^  '  '     ^     2    J-j^f ^ijCk^IF^'-^ 


As  was  -tjointed  out   before,    these  results  are  well  knovm. 


-3^ 

9.      The  Caae   of  a  Hoa-Unlfom  Atmogphera 

If  the  ataosphere   !•  not  uniform  and  is   such  that  r  x*'(r)   a«t   no  station.iry 
point,    that  Is,   phyflically,    the   atmosphere  doee   not  hsTe  a   duct,    then   It   can  "be 
treated  siailarly   to   the   ca«e  of  a  iinlform  atmosphere.      We   shall   see   that    the  results 
of  the  preceding  sections  will  "be   only  slightly  modified, 

first,  we   find  an  approximate   formula  for  T  as  defined  in  (^3).      Since   a  will 
\)e  near  to  r^   andX        to  ak(a),    the  bracket   in  (53)  may  be  approximated  as   follows. 
Put 


Then 


(67) 


2ak   ^  ^  T 

/v  a^  f  (r-r  )  -V  2  -2^   (  £.  .  1)  =  21^  (  ^  -  1; 


(x  a;   -'  X        0  0 

o  0 

e  r^k(rj  -A^^^  . 
rk(r)  -X^/2-  5  (r-r  ) 


Here,  we  have  made  use  of  the  fact  that  sine. 

»  00 


approxiaately , 
How, 


k  aT-  2^2^^  K  f  C  ^  -  1)^/2  ^  ^   23/2  j  j,  ^  (    r.  _  ^j3/2/ 

o  o    I  r  or  '  -'     . 

Jr  °  • 


If  ve  pat,   as  before. 


"'0    [^*   ^(Ec^*)-"/^] 
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than 

and  the  eigenralue  equation  (^5)  "beco«i«8  the  following: 

(6«)        (rk^ar^/3(2/)V2zlj^(2/p/2/3|  ..i  ^   y  1  -  k^  /  k^2  -  ^   , 

where  >?    le   defined  by   the  last   equation, 
Hov 


t^    X^''^  =  ^MrJ  -  ak.  ♦^(r  -a) 


o       o 


V  Ta/lor's  Theoren,   but 

=  -    J"  r     (Ec  a)"2/3/v    -a    (T  (Ik  aF^/^ 


r     -  a  .     _      . 

O  O  0 


and  BO 

£k 


(69)  t  =  ak.  -<r(Ec  aW3         0 

a  0  c 

a 

In  most  cases  the  term  ak  (a)  will  be  negliglbla  coatoared  to  k  ,  ao  that 
^  a 

L  =  k  and  I  =  k  /k  .   Than 
7    a         a'  o 

(70)  t  -  ak  -  <r(k  a)^/3  ^ 

a      a 

If  this  fomnila  for   t   Is  cwpaxed  with  the  corresponding  fomnila  (63)  for   the  case 

of  the  uniform  atmosphere,   we  observe  that   the  only  chang«  is   the  substitution  of 

k     for  k  . 
a  0 

The  value  of  M  can  be  found  by   the  method  used  in  the  preceding  section,      '.fe 

hove 
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■ow,   froB  (69), 


XT  =  -  E^  (^.-r^' 


10  that 


Again,   consider  the  ease  where   the  transmitter  and  receirer   are  both  on  the 
larth'e   surface  so   that  h  =   r  =   a.      Using  (69)   and  (71}   in  (U-6),  we  haTe 

in)  ^ ^ 

erpClk  d)        V     2niXi 
TT  a  9  ,     /_    ^-l  1    ^'^l' 


a       0 


^are 


J^,-(tta)V3e 
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Hote  that  egaln.  If  ak  Ka)   is  negligible  to  k  ,  go  that  K  =  k  /k 

"  a'  o   • 

7^^  =  (k^a)  ^/3  0 

r^  (/^)  =  ?(?<:) 

aa  defindd  in  the  preceding  section  and 

expClk  d) 

This   formula  shows   that  the  Herte  potential   in  the  case   of  a  non-nnifora  atmosphere 
Can  be   ohtained  from   the  Hertz  potential   In  the   casea   of  a  constant  atmosphere  ty 
■laking   the  suhsti  tut  ions:      k     for  k     and  ^  for  ^  where  J?  is  obtained  fron  r?  hy   re- 

placia^  k^  by  k^. 


The  case  where   the   transmitter  and  receiver  are  not   on   the  earth's   surface 
can  be   derived  frc«   the  preceding  resiilts  by  multiplying  each   term  in  the   series 
for  r(^)  by   the  corresponding  height-gain  factors  for   the  transmitter  and  receiver. 
The  height-gain  factor  for  the  receiver  le 

by   (5^).      '"'e   shall   introduce  into  this   formula  0    ,    the  angular   distance  of  the 
receiver  froo  the  horizon,   and  then   (75^  will   take  a  simpler   form, 

from  the  definition  of  9   ,    and  since  r  -  a  is   small  compared  to  a,   we  hAve 

COS  d     -     ie/l  -     ^^=^ 
r         r  a 


and  sine* 

>2, 


cos  9 


,-i-e;/2    . 
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we  h&y 


then,  from  (67) 
2 


e^  = )/  2U-a;/fl 


/^r      a      I 


(ff- 


-^^^  a;  2K^ 


0 


(t-  0"«'[H-^<-.'"'']-] 


K^(^.)-2/5(>i!-2f) 


and  ve  find  that 


Since  for  r  =  a,  %     ■  0,  we  see  that  the  rig)it  hand  laeober  of  (75)  ■becomes 


(76) 


■?^^^^. 


-J./ 2 


2/ 


i  J 


Hi)^{(;f2,3^p/2/3| 


Hi)][(2n3/^/3] 


There   is  a  similar   formula  for   the  height   gain  fivctor  for   the   trananltter.      Instead 
of  X       in  (7b)  we  must  useXv  wi^'^e 


;j/^.(ak//3o^,  e^»y/2(w)/a     . 


It   shoold  he  pointed  out   that   r^stilts   similar  to   those  of   thi«   section, 
which  show  that   the   field  due  to  a  non-cmiform  atmosphere  could  he  obtained  from  the 
field  in  a  uniform  atmosphere,   have  "been  obtained  by  Bckersley  {lOj   and  Breaaer  [llj  , 
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How«V9r,    th«y  had   to  asstiaa  a  particular    for*   of  vprlatlon  of   the   dlalectrlc   constant 
to  irlt, 

e(r)  -  1  -^*  Pa^/r^ 

wharee(    and  ^  ar«  amall   lonstants.      The  reaults    of  this   section  we  Indflpendant   of 
the  form  of  variatloa  of  the  dielectric  coastant. 
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10.     Cooaectlon  ylth  the   Tlat  Earth"  Theory. 


The  preceding   sections  have   shovn  that  to  a  good  order  of  accTiracy   the   eig-sn- 
depend  on  the  "behavior  of  k* 
point,    that  is,    the  point  at  vhieh 


values   depend  on  the  "behavior  of  k  (r)  only   in  the  neighborhood  of  the   transition 


k^(r}-     % 


r 


It  was   also  shown  that  at  tke  transition  point  r     =   a^  ver^'  closely.      If  these  facts 
are  taken  into  account,    the  differential  equation   (3U)  may  "be   approximated  as 
follows : 


(77) 


where 


'^*l<^-<]  <")-  ° 


H^ 


r\hr) 


is   the  modified  index  of  refraction  euid  where 


k^  -     t  (t  ♦D/a^     . 
n  n     n       '  • 


Xquation  (77 )   ^^  exactly  the  eigenvalue  equation  that  la   obtained  in  the 
tiropoephere  theoiy  on  the  a^sujiption  of  a  flat   earth  ri^J.      If  the  V.K.B.   method 
is  used  to  solve  this    equation,   we   see  that  we  get   the   same   results   as  by  applying 
the  V.K.B,  method  to   the  more   exax:t   equation  (3^).      There  will   be  a  discrepsmcy 
between  the  two  methods   only  at  points  wh«re  r  differs  very  much  from  a.      This   indi- 
cates  that  the  height-gain  factors   on  the    "flat   earth"   theory  will   be   wrong  for 
large  heights   above   the  earth's   surface.     Another  possible  source  of  error  is   due  to 
the  fact  that   the  Legandre  functions   of  the   exact   theory  have  been  approximated  by 
exponentials.      This   approximation  will  be  bad  for   extremely  large   distances   between 
recelrer  and  transBitter.      Similar  conclusions   on  the  range  of  validity  of    the   "flat 
earth"   theory  were   reached  bjT  Pak«rle  [7 J  .      On  the  whole,   however,    the  results   of 
the  prerioua   sections  have  shown  that    the   "flat   earth"   theory  gives   a  very  satis- 
factory approximation  to  the  solution. 


11,     An  Atmosphere  Containiag  a  Pact. 

la  Section  9  ^^  haT«  dlscueaed  propa^^ntlon  In  a  non-uniform  atmosphere  under 

2  2 

the  asiuaptlon  that  r  k  (r)  does  not  hnre  a  etatlonary  point  near  r«  a.   The  reason 

for  this  nsminptlon  Is  that,  otherv-lse,  the  differential  equation  (3U)  will  her*   two 
transition  points,  that  is,  the  equation 

k^(r)-  ^  -  0 

r 

will  have  two  roots   In  ^  near  X-    (alc^)   .      How,    the  Langer  approxLmntion  by  mesne   of 
Beaeel   functions   of   order  one-third  is  Yalid  in   the  neighhorhood  of   ench   transition 
point   separately  but   it  cannot   be  used  In  a  neighborhood  of  both  transition  points. 
Lander  \llj  haa  pointed  out   that   in  such  a  case   the  approiimRtion  to  the   solution* 
of  the   differential   equation  should  be  in  terns    of  confluent   hypergeonetr ic   functions, 
Since   these  functions   are  not   too  well  tabulated,   we  shall  attempt   another  approach. 
In  prectlce,    the   two   transition  points   of   the  diff ar(»ntlal   eouation  ore  close 
to   each  other.      We   ehsdl  approximate   the  differential  equation  by  nsBumln^   thnt  the 
transition  points   coincide  so  that  we  have  a  double   transition  point.      In  such  a 
case  we  have  ^  =   2   in  equation  (5I)   and  the  approximate    solution  will    be    expressed 
in  terms   of  Bessel   functions   of  order  one-fourtk. 


Pat 


Let  r     be  a  zero  of 

k^(r)  - 

^.0. 

k     a 
0 

;l 

(^%) 

k^  a2  r^ 
0 

0 

[ 


k^Cr).     4 


T"  r 
0 


and 


g(r)  -  (iN-r^)  I  k^r) 


r"  r 

o 


•0   that  f(r)  has   a  double   «ero  at   r  ■  r  .     According  to   (50), 


2dr 


and  then  froai  (52)  since  o(  =  2,  we  hure  that 
The  eigeoTalue  equation  nov  becoaee 


2 


In.  case  the  earth  Is  assuned  to  he  a  perfect  conductor,  the  equation  reduces  to 

H^/1^  (k.ar  )  =  0 


0  a 


The  first  two  roots  of  the  equation 

,(1) 


Bi)i   C7^)  =  o 


are  7^-   (2.38  •»■  O.I71)   e""^    ,    ^    «  (5.6  ♦  O.I71)  •""^  ,  so   that  the  eigen- 

ralue  equation  becomes 

(SO)  k:  a?   -   (2.38  +  O.I71)   ^""^    • 

0      m 

The  integral   for  7~  can  he   easiljr  approximated  in  certain  cases.      If 
f       (r   ) 


o 


then  we  may  write 


f(r)^|     f"   (r   )   (r-r   )^    . 
2  00 


lov 


f' 

\^ 

- 

i 

i,;i 

r« 

r 

0 

2? 

0 

2     U 
a     r 

J       ■ 

ty 

'  the 

last    equation. 

80 

thftt 

f(r)  « 

r 

(r-r^)2 

and 


J    T         ^  0 


Put,    as  "before 


a  »  r     (1  •►y  ) 

0 


in  equation  (81),   then 

and  in  case  the   earth  is  a  perfect  conductor,   we  Bust   solve   the  equation: 

I  (ra)V2j2.   (k^a)-^   (2.38  ^  O.iyi)  e"^      . 

Sote  that  ^  will   depend  on  r     and  therefore  also  on  ^  .      In  case^    '      Is   a  slowly 
Tarying  function  of  J*    ,   we  may  consider  it   a  constant  and  thea 

Since 

t    ^^^/^-k(r   )^'^'k     r-vk     a       (l-f) 

O  O  00  00 
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ve  find  that 


(82)  t   ^k^a-^   (k^/|  )^/2  (k^a;^/Ni+.76*  0.3^1)^/2     e""^/^  . 

ronaula  (82)    shows   that    the   ijaaglnary  part    of   t  varies    as      (k  a)   '       for   the 

o 

cftse   of  a  duct,   as   contrasted  to   the   case   of  a  hoaogeneous   atmosphere  where   It  T^rled 

as   (k  a)  '-^  ,      Since   the  attenuation  of  the   field   in  decibels   depends   on  the   Imaginary 
part   of  t,  we   can  see   that  the  range  vrill  he  greatly   Increased  in  the  case   of  a  duct, 
for  (k  a)  of   the  order   of  10   ,    the   range   is  multiplied  toy  h,6. 
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12,      Th«  Eectrlc  Dlpole. 

The  preceding  work  has  been  carried  out   for  the  case  of  the  magnetic   dipole 
tecanse   then  the   equations  and  the   discontinuity  condltione  beoone  fliLipler,      How- 
ever,   the  case   of  most  practical    interest    is   that    of  an   electric   dipole.      In   this 
sectioa  we  shall  give  hriefly  the   equations  appropriate  for  an  electric   dipole, 

Jirst,   a  new  Hertz  rector  and  potential  iraBt  "be  defined.      Instead  of   (2)  we 
put 


(83)  H»    VXC* 

+ 
where  C^^  again  has   only  its  radial   component   not   lero  so   that 

Scl  =   (rU^.   0.  0)     . 

Instend  of  (6),   we   find 


r  sin  e 


and  instead  of  (8),   we  have 


) 


^-^  =  -   ^   •    *e  =   -  IT-    973^   ^^1^  .    «r  -  =e  =  ^r  °   . 


Since  I-  and  H^  BTi-st  "be  continuous  across  a  discontinuity  surface  r  ■  constant. 

'^,1      1  3^1 
it  follows  that  U,  and  7-  _■  ■   naet  "be  continuous  across  this  surface.   Instead 

of  (lU)  we  have  the  following  equation: 


r   •^ 


with  the  same  boundary  conditions  as  for  (l^).   Put  a^in  t  (r)  for  those  solutions 

of  the  hoBogeneous  equation  corresponding  to  (86),  which  are  regular  at  the  origin, 

and  V   (r)  for  those  solutions  of  the  honogeneous  equation  which  satisfy  the  radiation 
n 

condition  at  infinity.  We  find  that,  now,  the  Vronskian  is  no  longer  a  constant, 
but  that 


-k&~ 


iSj)  ".^J^V  6-,(r)  -  constant. 

The  forvula  for  n.(r«0)  is  the  saae  a*  (I9): 


(r)w_(b) 


•  ^^'.^^  -  h  ZI    ^-n^ 


do. 


(88) 


1   * 


r  <  t 


r  <  "b 


Ih«  contour  integral  representation  and  all  the  transformations  are  unchanged, 
EoveTer,  the  eigenvalue  problem  becoses  the  following: 

yind  those  ralues  of  t  for  whioh  there  exists  a  solution  of  the  differential 
equation; 


(89) 


^hl'^   ^    U^.li-ii)(„,  =  o 


which  is  "both  raguleo"  at  r  =  0  and  satisfieB  the  radiation  condition  (I5)  at 
infinity. 

Because  of  the  different  discontinuity  conditions,  the  eigenvalue  equation  is 
changed.   Instead  of  (Ul)  we  have 


ftTTT   dF^^t^ 


rr. 


7r^{'^''V./a'vj 


r"  a 


'"'  \n/a<V^ 


r«  a 


80  that   the  eigenralue  equation   (W+)   hecoBes 

d 


(90) 


rw. 


r»  a 


.If.^ 


After   transforming  the  contour  integral  for   (88)   and  using  the  asyBptotic   formulas 
for   the  Legendre  functions,   we  get   the  following   result,    instead  of  (U6) 


-u^ 


r(2t,n)l-|^/2            1(^1/2)9  V       (b)  w       (r) 

(91)  U(r.O)  .  -    V-    -~f_—  _J2 1__ _J___!j^ 

fi«re  t    (j  ■  1,    2,    ...)   are   the  elgenvaLuee   of   (89), 


aquation  (89)   can  be  treated  by   the  W.I.B.   method  after  a   transformation  le 

.      Put 


■ada,    to  get   rid  of  the  6"     factor.      Put 


rr  »    €    '      rr 
then  eqTiation  (89)   becoaes 

dr 


(92)  l!^(rv)*(k2^     iiilli)    (r^)  =  o 


.2 


dx*^ 

This   transfomntlon  is  relid  only  if  we  aesune  6    is   continuous   and  differtntisble. 
Since  we   ahall   oixly  apply   the  V.K.B,   nethod  to   (89)   for  r  >  a,   our  assumption  is 
Justifiable, 

The   eigenvalue*   at   t.   of  (92)   caoi  be  found  In  exactly  the   same  way  as   we   found 
the   eigenraiues   for  (3^).      We  shall   not  repeat   the  work  but   one   difference   should 
be  noted.      In  case    the   earth  is   assumed  to   be   a  perfect   conductor,    eo    6  ■•    "*  ^    "    00, 
the  eigeoralue   equation  (9^)  becoaes 


»  0 
r«  a 


or 


i    «  r»  a 


instead  of 


rw       »  0     . 
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Appendlx  I  -  The  Contour  Iategr>-1  on  the  Semicircle  . 

In  thle  eectlon  we  conrplete  the  dlsrueelon  of  the  contour  Integral  represen- 
tation In  (29)  by  showing  that  the  Integrala  over  the  semi-circle  in  the  first  and 
second  quadrant  yanlshee  as  the  radius  of  the  circle  goes  to  infinity.   This  Justi- 
fies the  expansion  of  the  Integral  as  a  sun  of  residues  and  an  integral  over  the 
imaginary  aixis. 

Consider,  now,  the  Integral 


,   .        /"(t  +  l/2)y.(b)w  (r 

"        J ^~■^ 


)   p. (-cos  «) 

-    -^ — r-     d  t 

sin  TT  t 

10 


over  a  quarter  circle  [  1 1  «  T  in  the  first  quadrant.   Put  t  =  T  e   ,  then  from  (Uo) 
we  hare 


9T 

To  estimate   the  other  factor  in   the  Integrand,  we   shall  Investigate   the 
aeyraptotlc  behavior  for  large  t  of  the   solutions   of   the  equation   corresponding  to 
(16): 

jiSsX    .(,2.    iiiti))    („)     .0     . 
dr  r 

It  can  he  easily  verified  by  eubatitution  in  the  differential  equation  that  there  are 

two  Independent  solutions  which  behave  for  large  t  either  llk» 


V-n^-^/>^'] 


rv'^  C,r    II  -  :—r    I      Ic  r  1  dr 
or 

where  C,  and  C,  are  independent  of  r.   Since  rv  was  by  definition  regular  at  r  =  0, 

t+1 
it  must  behave  lllce  C.  r    for  large  t.   When  these  asymptotic  formulas  are  differ- 
entiated with  respect  to  t,  we  find  that  for  all  solutions  of  the  differential 
equat Ions 
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Ik -I    .. 


Vow 


(t+l/2)T^(l))w^(r)  (VH/2)bT^(\))rw^(r) 


t+1/2 

bT^(b) 

br 

M 

where  d   /        v  d   /        v 


K  - 


\  "t 


M  will  vanish  If  and  only  If  W.    =0.      If  T  i«   eo  chosen  that   the  q-oflrter  circle 
does  not  go   throtigh  einy  s«ro  of  W   ,    then  M  will  he  bounded  below  by  some  multiple 

t  A 

of  t/r.     We  have  then 

(t^l/2)v^(b)w^(r)  ct/b\^^       1  C/b\*-^^ 


^ 


Ct/b\  -^       1       ^     C/b\ 
br  Vr/  t/r  b  V  r  / 


and  since  r  <  b,  this  tera  will  go  to  zero  as  T  goes  to  infinity.  Tlnally,  the 
integral  will  be  less  in  absolnte  value  then 


M     ^2    m^^«       /-wxT  cos  0 


b«T 


and  as  this  clearly  approaches  zero  for  large  T,  the  shifting  of  the  contour  in 
section  \  was  Justified, 
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4pT)«nAix  II.  An  Estimate  for  the  Integral. 

It  hae  been  ehown  that  the  contour  Integral  repregentatlon  of  the  Hertt  poten- 
tial can  he  expressed  (39)  as  a  eum  of  reBldues  plus  eux   integral  over  the  imaginary 
axis.   In  this  section  we  estinate  the  value  of  the  integral  and  ahow  two  things: 

Tirst,  the  integral  Taniehes  identically  if  the  earth  is  assumed  to  be  a  per- 
fect conductot;   second,  the  intefi:ral  is  small  compared  to  the  first  term  in  the  i 
residue  series  if  the  complex  dielectric  constant  of  the  earth  is  small  comi>ared  to 
the  dielectric  constant  of  the  atmosphere.   Since  this  condition  Is  always  satisfied 
in  practice,  we  are  Justified  in  neglecting  the  contribution  of  the  integral. 

In  order  to  estimate  the  Integrand  in  (39)  we  must  make  use  of  the  fact  that 
the  differential  equation  (3**)  !•  unchanged  when  t  is  replaced  by  -t-1  so  that 

»r^  (r)  -  w_2i  (r)      and      w^^^r)  -  ^[^^   (r)   .  • 

Since  w^  &ai.   w}   are  independent  solutions  of  (3'^)>  we  nay  write 

t  t  f 

T^(r)  .c<w^^^    (r)  +^wj2^    (r)  | 

! 

whereat   and  (3  are   independent  of   r  but  may  depend  on  t.  j 

Hote   that  | 

w^     -^Wronsilan   (w[^^      .  w^^^    )     =  ^ik^  ^    • 

2    2  i 

Now,  we  make  use  of  the  fact  that  for  r  <  a,  k  ■  k.  ,  a  constant,  so  that 


rT.(r)     -     r^^^  J*^Wo    (l^.r)  yr<a. 


t^^^  ^        "t+i/a  ^''l 

It  the  boTindary  r  >  a.  we  mist  hare 

.(1)    /.^     ^a,j2)    f.s  .1/2 


O^awJ^'    (a)     +^aw;^'    (a)  -  a^^'  Vl/2^V^ 

<h^<%   ^^h^<\'   h^^'^'^^Uz^l 
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Th«8e   two  eqviAtlont  will  d«t«rmine  a(  and  A  ,     We  find 

rw^2)  d^/    1/2   .  \         1/2-  d        ,      (2)    V 

^     .  ^.        ^^     .       '^t-Kl/2)-  ^       Vl/2     "d?  ^"^t        ^     ■ 

■  -  — rr\ **^  t   >^eB  M.  ie  defined  "by  tL»  last  equation, 

aw^   (a) 

Nov,  consider  the  expression 

.  ''''"■>->'"'-'        •^''<-'      (!)<,,  (1),  ,  « 

o  t 

lyerythln^  in  this  formula,  except  the  factor  M.  ,  is  tmchanged  when  t  is  replaced  "by 
~t-l,  80  that 

-t^^K^'-)      ^.t.l^^>  -.t-l^')        wf  ^a)wjl^b)w^^^r) 


How 


(S6) 


W  V 

t  "-t-1  2i)c  w 

0 


V 


♦  1 


and  we  aee  that   In  thie  formla  the   only  term  whldi  changes  wheo   t  la   replaced  hy 
-t-1    19 


t?Ll<«(-'''Vl/2>]. 


Hotrerer,    the  argruaent   of   this  Bessel  function   is  k,   a,   which  1b  a  rery  large   t^-aan- 
tity  since  the  conductlrlty  of   the  earth   is   so   large,   and  hy  the  asyatptotic  formula 
for  the  Bessel  functions  of  large  argtunent  we  hare 

Hote  that  this  formula  Is  TStlld  only  If  k.a  Is  Ysry  much  larger  than   t.   How  aince 
k.    Is   larger  than  k,    the  fraction  In  the  formula^niuBt  be  very  small  so  that 

H.    -  ^^fi.]   i«  nearly  zero. 

I 
If  we  put   t  equal   to  lT-l/2  so  -t-1  =  -cT'    •  l/2,   we  see   then  that  the 
■bracket  in  the  Integrand  of    (39)  nreuit  he  nearly  sero  for  T'very  much  smaller  than 
k.,a.      In  the   remaining   interra^  of   integration  we  use   the   formula   (^0).     We  find  then  that' 
that   the    integrajid  Is   of  the  order  of  aoagnitude  e  .     Vow 


/ 


00 


-ro  ^«.      1     -k»« 


1/3  ^  l/v- 

and  we  see  that  this  is  rery  small  compared  to  e  or  to  a    '      which 

is  what  Is  obtained  from  the  residue  series.  We  have  thus  shown  that  the  Integral 

In  (39)  is  small  compared  to  the  first  term  in  the  residue  series. 
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